Radio frequency spectroscopy and the pairing gap in trapped Fermi gases 
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We present a theoretical interpretation of radio-frequency (RF) pairing gap experiments in trapped atomic 
Fermi gases, over the entire range of the BCS-BEC crossover, for temperatures above and below T c . Our 
calculated RF excitation spectra, as well as the density profiles on which they are based, are in semi-quantitative 
agreement with experiment. We provide a detailed analysis of the physical origin of the two different peak 
features seen in RF spectra, one associated with nearly free atoms at the edge of the trap, and the other with 
(quasi-)bound fermion pairs. 



PACS numbers: 03.75.Hh, 03.75.Ss, 74.20.-z 

A substantial body of experimental evidence for superflu- 
idity in trapped fermionic gases (jj, |2j, yl |4[ has focused at- 
tention on an important generalization of BCS theory associ- 
ated with arbitrarily tunable interaction strengths; this is called 
"BCS-Bose-Einstein condensation (BEC) crossover theory" 
|3l • This tunability is accomplished via magnetic field sensi- 
tive Feshbach resonances. At weak interaction strength con- 
ventional BCS theory applies so that pairs form and condense 
at the same temperature, T c , whereas as the attraction be- 
comes strong, pairs form at one temperature (T*) and Bose 
condense at another (T c < T*). The intermediate or uni- 
tary scattering regime, (where the fermionic two-body s-wave 
scattering length a is large), is of greatest interest because it 
represents a novel form of fermionic superfluidity. In contrast 
to the BEC case, there is an underlying Fermi surface (in the 
sense that the fermions have positive chemical potential fi), 
but "pre-formed pairs" are already present at the onset of their 
condensation. 

The difficulty of obtaining phase sensitive probes and the 
general interest in this novel superfluidity make experiments 
which probe the fermionic excitation gap extremely impor- 
tant. While in the weak coupling BCS limit the gap onset 
appears at T c , in the unitary regime this gap (or "pseudogap") 
appears at a high temperature T* and directly reflects the for- 
mation of (quasi-)bound fermion pairs 01 01- For the 
trapped Fermi gases, one has to devise an entire new class of 
experiments to measure this pairing gap; traditional experi- 
ments, such as superconductor-normal metal (SN) tunneling 
are neither feasible nor appropriate. The first such experiment 
was based on radio frequency (RF) spectroscopy |7]; this fol- 
lowed an earlier proposal by Kinnunen et al |9, 10], who also 
presented an interpretation of recent data in 6 Li II ill in the 
unitary regime. However, some issues have been raised about 
their interpretation in the literature fl2ll . Moreover, the spectra 
in the BEC and BCS regimes also need to be addressed. 

It is the purpose of the present paper to present a more sys- 
tematic analysis of RF pairing gap experiments for the en- 
tire experimentally accessible crossover regime from BCS to 
BEC, as well as address recent concerns fl2ll . Our studies 
address the two peaks in the spectra observed experimentally 
at all temperatures, and clarify in detail their physical origin. 
Essential to the present approach is that our calculations are 
based on trap profiles 11311 and related thermodynamics 
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which are in quantitative agreement with experiment at uni- 
tarity Il4ll5ll where there is a good calibration. 

In the RF experiments 0], one focuses on three different 
atomic hyperfine states of the 6 Li atom. The two lowest states, 
|1) and |2), participate in the superfluid pairing. The higher 
state, |3), is effectively a free atom excitation level; it is un- 
occupied initially. An RF laser field, at sufficiently large fre- 
quency, will drive atoms from state |2) to |3). 

As in Refs. we base our analysis on the conven- 

tional BCS-Leggett ground state 0, extended |5, 6] to ad- 
dress finite temperature effects and to include the trap poten- 
tial. In this approach, pseudogap effects are naturally incorpo- 
rated. We begin with the usual two-channel grand canonical 
Hamiltonian H — /j,N 11811 which describes states |1) and |2), 
as in Ref. 1 16], and solve for the spatial profiles of relevant 
physical quantities. As a result of the relatively wide Feshbach 
resonance in 6 Li, the fraction of closed-channel molecules is 
very small for currently accessible fields. Therefore, we may 
neglect their contribution to the RF current, as was done in 
Ref. H. 

The Hamiltonian describing state |3) is given by Hj, — 

M3^3 = Lk( £ k + ^23 - M3)c3,k C 3,k' where e k is the 

atomic kinetic energy, c 3 k is the annihilation operator for 
state |3), ^23 is the energy splitting between |3) and |2), 
and fj,^ is the chemical potential of |3). In addition, there 
is a transfer matrix element Tk iP from |2) to |3) given by 

Ht = Ek P (^ k .p c 3.p c 2k + ^- c -) F° r pl ane wave states, 
Tk, P = T5(q L + k - p)<S(w kp -Wi). Here q L « and u> L 
are the momentum and energy of the RF laser field, and oj^ p 
is the energy difference between the initial and final state. It 
should be stressed that unlike conventional SN tunneling, here 
one requires not only conservation of energy but also conser- 
vation of momentum. 

The RF current is defined as I = {N 3 } = i{[H,N 3 ]). 
Using standard linear response theory one finds / = 
2T 2 lm[X ret (— ojl + fJ-3 — m)]- Here the retarded re- 
sponse function X ret (uj) = X(iuj n — ► oj + i0 + ), 
and the linear response kernel X can be expressed in 
terms of single particle Green's functions as X(iu> n ) — 
T E m ,k G 3(k, iv m )G(k + q L> iv m +iu n ), where w„ and v m 
are even and odd Matsubara frequencies, respectively. (We 
use the convention h = ks = 1). After Matsubara summa- 
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tion we obtain 

=i2 



BEC 



Unitary 



BCS 



1 = T- \ &v E A 3(k, v)A{ v , i/)<5(q i+ k-p) [/(„') - /(„)] 

^ J k,p 

(1) 

where i/ = v — ujl+hz — h, and f(x)is the Fermi distribution 
function, ^(k, j/) = 27T(5(^ — (e k +W23 — ^3)) and A(k, !/) = 
— 2ImG(k, ^ + i0 + ) are the spectral functions for state |3) 
and j 2), respectively. Finally, one obtains 



I(u) 



2tt 



^2 A(k + q L , e k - u - (j,) 

k 

x [/(e k -u-fi)- /(e k + w 23 - fJ.3)] , (2) 



where = ujl — U23 is defined to be the RF detuning. 

We have introduced a T-matrix formalism for address- 
ing the effects of finite temperature based on the standard 
mean field ground state @, @]. Here the Green's function 
G(k, v) contains two self-energy effects deriving from con- 
densed Cooper pairs as well as from finite momentum pairs, 
(which are related to pseudogap effects). These finite lifetime 

A 2 

pairs have self energy T, pg (k, v) w v+€ ^ l+ , n where 7 ^ 0. 
By contrast, the condensate which depends on the superfluid 

A 2 

order parameter (OP), A sc , enters with S sc (k, v) — ly+e °l M - 
The resulting spectral function, which can readily be com- 
puted from £ = Y, pg + S sc , is given by 



A(k, v) 



(^ + a)> 2 



e 2 

Sk 



Here £ k = e k - \x. E k = 



(3) 



A 2 (T) is the quasiparticle 
dispersion, where A 2 (T) = A 2 SC (T) + A 2 g (T). The pre- 
cise value of 7, and even its T-dependence is not particularly 
important, as long as it is non-zero at finite T. As is consis- 
tent with the standard ground state constraints, A pg vanishes 
at T = 0, where all pairs are condensed. It is reasonable 
to assume that 7 is a monotonically decreasing function from 
above T c to T = 0. Above T c , Eq. (0 can be used with 
A sc = 0. Because the energy level difference W23 (~ 80 
MHz) is so large compared to other energy scales in the prob- 
lem, the state |3) is initially empty. It is reasonable to set 
/(e k + ^23 - M3) = in Eq. (0. 

For the atomic gas in a trap, we assume a spherically sym- 
metrical harmonic oscillator potential V(r) — mu) 2 r 2 /2 in 
our calculations, where u> is the trap frequency. The density, 
excitation gap and chemical potential will vary along the ra- 
dius. These quantities can be self-consistently determined us- 
ing the local density approximation (LDA). Here one replaces 
/1 by a spatially varying chemical potential /i(r) = [i — V(r). 
The same substitution must be made for /13 as well. At each 
point, one calculates the superfluid order parameter A sc (r), 
the pseudogap A pg (r) and particle density n(r) just as for 
a locally homogeneous system; an integration over r is per- 
formed to enforce the total particle number constraint. Equa- 
tions (|2ji and l|3} can then be used to compute the local cur- 
rent density J(r, u>) and then to obtain the total net current 
I{uj) = J d 3 rl(r,uj). 
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Figure 1: RF spectra I(uj) for the near-BEC (720 G, left column), 
unitary (837 G, fc^a = oo, middle column), and near-BCS (875 
G, right column) cases as a function of RF detuning uj. The values 
of T (except in the first row) and ifo were chosen to match the 
experimental values in Ref. 



Figure^shows the calculated RF excitation spectra in a har- 
monic trap for the near-BEC (720 G), unitary (837 G, a = oo) 
and near-BCS (875 G) cases, from left to right, for a range 
of T from above to below T c . Here, for definiteness, we take 
•y/Ep = 0.1 + 0.3T/T C , which monotonically increases with 
T as one may expect 1 19]. The values of T used in all rows but 
the first (which involved the less interesting extreme Boltz- 
mann regime) were chosen to be consistent with the corre- 
sponding values of T" used in Ref. |7], on the basis of a the- 
oretical thermometry 11611 . Here T' refers to the initial tem- 
perature of an isentropic sweep starting from the BEC side 
of resonance. The values of kpa were calculated from the 
known values of Tp and a. Just as in experiment 01, two 
distinct maxima are seen. A very sharp peak at u> = ap- 
pears only for T ^ 0; this peak is, thus, related to thermally 
excited fermion states. A second and broader maximum is 
present at sufficiently low T and is connected to the breaking 
of fermion pairs between states |1) and |2), with subsequent 
transfer of state |2) to |3). The broadening of zero uj peak, as 
T decreases, reflects the increasing values of the gap A. The 
general features of the spectra are in reasonable agreement 
with experiment for all three cases shown. 

The near-BEC plot is still far from the true BEC limit where 
kpa is arbitrarily small. Nevertheless, one can see from the 
lowest T figure that the absorption onset is only slightly larger 
(~ 5Ep as compared to ~ A.QEp) than the estimated two- 
body binding energy ti 2 /ma 2 , as expected. This near-BEC 
figure makes it clear that pairing effects are absent at the high- 
est T = l.OTp (> T*), where the free atom peak is sym- 
metric and there is no sign of a shoulder; this case is close 
to unitarity largely because of the size of kp. It is also clear 
from the middle figure that the "pairing gap" forms above T c , 



3 




*- — — ♦ * 



-20 20 40 

RF detuning (kHz) 

Figure 2: Comparison of calculated RF spectra (solid curve, T c ps 
0.29) with experiment (symbols) in a harmonic trap calculated at 
822 G for the two lower temperatures. The temperatures were chosen 
based on Ref. [7]. The particle numbers were reduced by a factor of 
2. 
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Figure 3: Origin of the two RF spectral peaks. Spatial distribution 
of (a) A and density n(r), contributions to (b) the low frequency 
peak at lo/Ef = —0.01 and (c) the high frequency peak at uj/Ef = 
0.33 at unitarity in a harmonic trap, calculated at T — 0.237V . (d) 
Schematic transitions from state |2) to |3), showing the extended k- 
space contributions to the low energy peak, at the trap edge where 
A < T. Here Rtf is the Thomas-Fermi radius for a noninteracting 
Fermi gas. 



as is expected. Although not shown here, we find that for the 
unitary case, there is an analogous pseudogap effect which ap- 
pears above T c via a shoulder in the spectra to the right of the 
oj = peak. Only when T > T* will this shoulder entirely 
disappear. At unitarity, we find T* w 2T C . Additionally, it 
should be stressed that the near-BCS case is still very far from 
the weak coupling BCS limit. 

Experimentally 1 20], one defines the (averaged) "pairing 
gap", Arf, as the energy splitting between the maximum 
in the broad RF feature and the to = point. For the near- 
unitary case in 6 Li (at 822 G), A RF /E F w 0.25 at the in- 
termediate T' = Q.5T F , whereas at the lowest T this ratio 
is around 0.35. The ratios found theoretically are roughly 
0.35 and 0.38 for these two cases. However, when the field 
is increased to precise unitarity (837 G) the numbers appear 
to considerably smaller with a ratio of < 0.2. On general 
grounds one can argue that very little change is expected with 
these small changes in field near unitarity. Anharmonicity as- 
sociated with a shallow Gaussian trap may explain this small 
discrepancy, along with possible uncertainties in the particle 
number ll2ll I22I1 . There may also be some interference with 
the Feshbach resonances between states |1) and |3) and be- 
tween 2} and |3) 12311 . which overlap with the resonance be- 
tween |1) and 1 2) but are not included in the theory. 

In Fig. [2] we compare our calculated spectra near unitar- 
ity (solid curve) with experiment (symbols) at 822 G for the 
two lower temperatures. The dashed curve is a fit to the data, 
serving as a guide to the eye. As in Fig.Q] we calculated kpa 
using the experimental values of Tp and a 112011 . After reduc- 
ing the particle numbers by a factor of 2, as in Ref. ll22ll . this 
brings the theory into very good agreement with experiment. 

To fully understand the RF excitation spectra, it is impor- 
tant to determine where, within the trap, the two frequency 



peak features originate. Figure |3ja) shows a plot of the den- 
sity profile n{r) and excitation gap A(r) at unitarity and 
T « 0.23T F « 0.83T C as a function of radius. Figures^) 
and|3jc) indicate the radial dependence of the local current 
I(r, u>) for (b) frequency near zero, where there is a sharp 
peak, and for (c) frequency near the pairing gap energy scale, 
where there is a broad peak. Just as conjectured in previous 
papers 171 ll ill , it can be seen that the low frequency peak is 
associated with atoms at the edge of the trap. These are es- 
sentially "free" atoms which have very small excitation gap 
values, so that they are most readily excited thermally. By 
contrast, the pairing gap peak is associated with atoms some- 
where in the middle of the trap. 

One might expect a rather broad free atom peak, reflect- 
ing a range of values of A(r) at trap edge, but this peak is, 
in fact, quite sharp, as is its experimental counterpart. The 
sharpness of the free atom peak is addressed via the schematic 
diagram of Fig. |3jd). When A < T (as in the trap edge re- 
gion), the dispersion of state \2) reduces to a simple parabola 
as for free fermions; it is thus similar to that of state |3), as 
seen in Fig. |3jd). Momentum conservation leads to vertical 
transitions shown by the arrows on the figure. It is important 
to contrast this picture with the situation for SN tunneling; 
here Pauli blocking effects are absent since the final state is 
empty for all k. As a result there is an extended volume of 
/c-space contributing to the transition at ujl = UJ23 (which 
corresponds to detuning u> = 0), thereby leading to the sharp 
spectral peak. At the very high T which were probed experi- 
mentally in Ref. [7], the sharp 10 = peak results in a similar 
fashion, although in this Boltzmann regime, the gap A is com- 
pletely irrelevant. 

A plot analogous to Figs. |3jb) and|3jc) can be made for the 
near-BEC case as well, to determine where in the trap the RF 
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gap Arf arises. It is easy to see that the threshold region is as- 
sociated with the trap edge where A(r) is small. Indeed, when 
[i < 0, the excitation gap is given by n 2 + A 2 . This implies 
that the two-body binding energy Ef, w — 2/i + airnh 2 jm sets 
the scale for this threshold, in much the same way as found 
for the deep BEC where a — > 0. As long as fi < 0, the 
values of A j^f and Ej, are very close, becoming equal when 
A — > at the trap edge. This supports the more detailed two- 
body analysis of these threshold effects in the BEC presented 
in Ref. 12411 . However, it should be stressed that there is an 
intrinsic rounding around the threshold, as seen in the bottom 
left panel of Fig. [2 

At T = these LDA-based RF calculations can be com- 
pared with the results lll2ll of Bogoliubov-de Gennes (BdG) 
theory. It should be noted that BdG theory is appropriate for 
the particular mean field ground state under consideration, but 
it cannot be applied at T ^ 0, since it does not take into ac- 
count the noncondensed bosonic degrees of freedom. A com- 
parison presented in Ref. 1 12] between a BdG calculation and 
its LDA approximation showed a difference in the low fre- 
quency tunneling current at T — in the fermionic regime 
(/j > 0). The finite spectral weight at precisely co = in the 
BdG result was interpreted to arise from Andreev bound states 
10. 

It was also speculated that at finite T, Andreev effects may 
be playing a role so that the free atom peak is possibly of a 
different origin from that considered here and elsewhere |7, 
HHl . It was noted in Ref. 1 12] that the BdG equations show 
that the entire trapped gas is in the superfluid state below T c , 
with A sc (r) being finite everywhere n(r) is finite. Therefore, 
it was presumed that the free atom peak found in Ref. full was 
an artifact of the LDA at T ^ 0, since in this approximation, 
there is a region of the trap where A sc = 0. 

In support of the present viewpoint it is important to note 
that the free atom peak derives from states where A(r) < T. 
The gap A, is the important energy scale, not the order param- 
eter A sc , for characterizing fermionic single-particle excita- 



tions. This can be seen from the fact that the spectral function 
of Eq. which enters into the RF calculations depends on 
A through E^, and is not particularly sensitive to A sc . From 
Fig. 0a) it follows that A is finite wherever n(r) ^ 0. It be- 
haves similarly to the nonvanishing order parameter in BdG- 
based calculations. Furthermore, both the OP in Ref. 1H2I1 and 
A (for the present case) behave as A ~ exp(— n/2kp\a\), be- 
coming exponentially small at the trap edge. Thus, as a result 
of pseudogap effects (which serve to distinguish A and Ay at 
any finite T) we believe that the concerns raised earlier 1U2II 
about the applicability of LDA for addressing RF experiments 
are not warranted. 

The results of this paper support a previous theoretical in- 
terpretation 1 1 1 ] of RF experiments [7] in the unitary regime, 
which applied the pseudogap based formalism of the present 
paper, albeit with approximated spatial density and gap pro- 
files. The present calculations avoid these approximations, 
and lead to spatial density profiles 1 13] and related thermody- 
namics which are in good quantitative agreement with 
experiment [13]. Our work clarifies the origin of the two 
generic peak structures seen in RF experiments, and addresses 
the entire magnetic field range which has been studied exper- 
imentally. The zero frequency peak comes from atoms at the 
edge of the trap, where A < T. We interpret these as "small- 
gap" rather than "in-gap" excitations. We have shown that the 
sharpness of this peak is associated with an extended momen- 
tum space available for the u> = excitations. The broader 
peak derives from the breaking of pairs and, except in the 
extreme BCS limit, this peak is present above T c , reflecting 
pseudogap effects. 
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